Angle dependence of Andreev scattering at semiconductor— superconductor interfaces 
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We study the angle dependence of the Andreev scattering at a semiconductor-superconductor 
interface, generalizing the one-dimensional theory of Blonder, Tinkham and Klapwijk. An increase of 
the momentum parallel to the interface leads to suppression of the probability of Andreev reflection 
and increase of the probability of normal reflection. We show that in the presence of a Fermi 
velocity mismatch between the semiconductor and the superconductor the angles of incidence and 
transmission are related according to the well-known Snell's law in optics. As a consequence there is a 
critical angle of incidence above which only normal reflection exists. For two and three-dimensional 
interfaces a lower excess current compared to ballistic transport with perpendicular incidence is 
found. Thus, the one-dimensional BTK model overestimates the barrier strength for two and three- 
dimensional interfaces. 
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I. INTRODUCTION 

An electron-like quasiparticle incident on a normal 
conductor-superconductor (NS) interface from the nor- 
mal side may become Andreev reflected into a hole-like 
quasiparticle with reversal of the signs of all three velocity 
components (retroflection) and of the energy (relative to 
the Fermi level) as shown by Andreev [jjj. Later, Blonder, 
Tinkham and Klapwijk (BTK) Q calculated the scatter- 
ing probabilities at a NS interface within a model where 
the scattering at the interface was represented by a delta- 
function potential barrier. The calculations were based 
on the Bogoliubov-de Gennes (BdG) formalism ||, for a 
one-dimensional (ID) geometry thus ignoring all effects 
due to quasiparticles with a momentum parallel to the 
interface. 

The BTK model has been widely used by experimen- 
talists to model normal-metal-superconductor junctions, 
and it has despite of its inherent approximations been 
quite successful in describing the main features of these 
devices. The quality of the junction interface has con- 
veniently been parametrized in terms of the normalized 
delta-function barrier strength. 

A more complete theory was developed by Arnold 
jlj using non-equilibrium Green's function techniques. 
The theory by Arnold furthermore takes the three- 
dimensional (3D) nature of the interface into account. 
However, the resulting expressions are complicated and 
require substantial numerical work. Generalization of the 
BTK model to tunnel barriers other than delta-function 
scattering potentials has been done by Kupka in a num- 
ber of papers P|§. Recently, Kupka M generalized the 



more realistic tunnel barrier model to include the angel 
dependence of the scattering. He found that by treating 
the scattering problem in the correct three dimensional 
picture, the effective Andreev scattering is reduced and 
the normal scattering probability is enhanced. Chaud- 
huri and Bagwell Q] and De Raedt, Michiclscn and Klap- 
wijk pj have also considered the angle dependence in 
their applications of the BdG formalism to the transport 
properties of NS interfaces. However, except for the ID 
work of Blonder and Tinkham [^| the above mentioned 
papers all focused on the case where there is no mis- 
match between densities (and hence Fermi wavelengths) 
or between effective band masses of the two materials 
forming the NS junction. In the case of SNS junctions, 
Kupriyanov (To) included effects of the parallel degree 
of freedom and different Fermi velocities of the N and 
S regions in his application of the Eilenberger equations 
to the dc Josephson current in junctions with clean in- 
terfaces. Using the BdG formalism, the effect of differ- 
ent Fermi velocities and effective masses was also consid- 
ered by Schiissler and Kummel jll] and Chrestin, Mat- 
suyama and Merkt jl2| in their numerical studies of the 
dc Josephson current in Nb-InAs-Nb junctions. 

Since much of the development in the recent years has 
been in structures where superconductors are combined 
with semiconductors, the goal of this paper is an analyti- 
cal study of the importance of the different quasiparticle 
propagation in the two materials, when the degrees of 
freedom parallel to the interface and effects of the un- 
equal Fermi velocities and Fermi wavelengths are taken 
into account. This is motivated by the observation that 
Andreev scattering cannot occur above a critical angle 
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where the momentum can no longer be conserved. The 
critical angle depends on the ratio of the carrier density 
of the semiconductor to the density of the superconduc- 
tor. Therefore one may expect larger differences between 
ID and 2D or 3D junctions for the case of a finite Fermi 
wave vector mismatch, which is indeed what we find. 

The effect of the angle dependence of the Andreev scat- 
tering probability is however somewhat suppressed by the 
fact that the current is carried mostly by particles in- 
coming perpendicular to the interface in 2D or at angle 
of 45° in 3D. Therefore we suggest an experiment where 
the angle dependence of Andreev scattering is probed in 
a more direct fashion, namely a mesoscopic device which 
explores the ballistic motion of quasiparticles and where 
the angle of incidence can be varied. Such a device is 
possible due to the advances in fabrication of mesoscopic 
semiconductor-superconductor interfaces (see e.g. Q), 
which have made it possible to study Andreev scattering 
in the ballistic regime. 

The paper is organized as follows: In Sec. II the BdG 
formalism is introduced, and in Sec. Ill the scattering 
probabilities at the interface is calculated. These scatter- 
ing probabilities are used in Sec. IV to calculate current- 
voltage characteristics and related quantities. In Sec. V 
an experiment is suggested and, finally, in Sec. VI dis- 
cussions and conclusions are given. 



II. THE BOGOLIUBOV-DE GENNES 
FORMALISM 



The BdG equations 



Wo(r) 
A*(r) 



A(r) 

-n* (r) 



(1) 



provide a microscopic formalism for studying inhomogc- 
neous superconductors and NS interfaces ||. Here, A(r) 
is the superconducting order parameter and 7io (r) is the 
Hamiltonian. In a general non-equilibrium situation, the 
Hamiltonian includes either a time-dependent vector po- 
tential or a spatially dependent scalar potential. How- 
ever, we follow BTK [|j and neglect the effect of a finite 
bias on the scattering probabilities which is justified if 
the height of the tunnel barrier is much higher than the 
applied voltage and/or energy of the carriers (relative to 
the Fermi level) [||. For an interface where the position 
of the conduction band and the effective mass change 
across the interface, we use the effective mass approxi- 
mation 14 15|] 



Wo(r) 



2m*(r) 



V + U(r) - /i 



(2) 



where U(r) is total electrostatic potential, and \x is the 
chemical potential. This approximation describes the 
spatial dependence of the dispersion relation, and the 
form of the Hamiltonian ensures conservation of the 



probability-current. For a discussion of justifications of 
this approach, see Refs. Jl4|,[l5| and references therein. 
We assume a parabolic dispersion so that the effective 
mass to* does not depend on energy (or momentum). 

The solutions to Eq. p) are vectors in the so-called 
electron-hole space (Nambu space), ip(r) = (u(r), v(r)) , 
where u(r) is the electron-like quasiparticle amplitude 
satisfying an ordinary electron-like Schrodinger equation 
and v(r) is the hole-like quasiparticle amplitude satisfy- 
ing a time-reversed Schrodinger equation. In electron- 
hole space, a probability-current-density can be associ- 
ated with the wave function, and is given by PJj|,|l4[ 



(3) 



J„ = film <^ u*(r) — Vn(r) - v*(r) —Vv(r) 

p \ y J m*(r) y ' K 777* (r) v ; 



The BdG equations and the conservation of the 
probability-current-density form the basis for our treat- 
ment of scattering of quasiparticles at the NS interface. 
Eq. (1) is used in calculating scattering amplitudes and 
the corresponding scattering probabilities are found us- 
ingEq. (3). 



III. SCATTERING OF QUASIPARTICLES AT A 
NS INTERFACE 

We consider a planar NS interface lying in the icy-plane 
at z = with a semi-infinite non-superconducting ma- 
terial for z < and a semi-infinite superconductor for 
z > 0. The superconducting order parameter is assumed 
to vary in space only along the z-direction. In order to 
solve the BdG equations, we include only scattering at 
the NS interface. Following BTK, we model the scatter- 
ing at the interface by a delta-function potential 



U(t) = H6(z) 



(4) 



where H is the strength of the potential barrier. For 
simplicity we neglect the phase of the pairing potential 
since only the absolute value is important for the con- 
sidered geometry. Furthermore, to avoid self-consistent 
calculations, we take the superconducting order param- 
eter to be zero in the normal conductor and uniform in 
the superconductor, i.e. 



A(r) = A 6(z), 



(5) 



where Ao is the BCS value of the energy gap and Q(z) 
is a Heaviside function. Similarly for the effective masses 
of the two materials, we assume that the mass changes 
abruptly across the interface 



l*(r) = 777 (N) e(-z)+ 777 (S) e(z) 



(6) 



where 



and 



are the effective masses of the 



normal conductor and the superconductor, respectively. 
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Eqs. represent the simplest forms of U(r), A(r) 

and m*(r) still capturing the main physics of the NS in- 
terface. 

Due to the simple form of the NS barrier potential, 
the superconducting order parameter and the effective 
mass, we can separate the variables and express the 



solutions in the parallel direction as plane waves, i.e. 
V> <N ' s) (r) = exp [i (k ( x N - s) x + k^ s) y)] ^ <N ' s) (z), where the 
superscript (N, S) refers to the non-superconducting or 
superconducting sides, respectively. Substituting this 
Ansatz into Eq. (|l|), yields the effective BdG equations 
for the z-direction 



2m< N »S) dz 



; + hs(z) ^ 
A e(z) 



ft 2 



2m( N . s > dz^ 



A e(z) 



« <N ' S) 



ip(z) = Etp(z) 



(7) 



where the effective chemical potential is defined as 

S) - m <n ' S) - ^ ([*r B, ] a + K s> ] 2 ) ■ (8) 

Eq. (Q) is mathematically identical to the ID BdG 
equations considered by BTK, and, therefore, we expect 
similar results for the scattering probabilities. We adopt 
the notation of BTK @, app. A] and all formulae for the 
eigenstates, scattering states, wave vectors (q ,± and A^) 
etc. are equivalent to those of BTK, but with the im- 
portant difference that the chemical potential is replaced 



by an effective chemical potential which depends on the 
parallel momentum according to Eq. (g) . 

We follow BTK and consider an electron-like quasi- 
particle incident on the NS interface from the normal 
side. At the interface it has an amplitude a of undergo- 
ing Andreev reflection, b of normal reflection, c of normal 
transmission and d of Andreev transmission. The scatter- 
ing amplitudes are obtained by matching the scattering 
states at the NS interface, using the appropriate bound- 
ary conditions for a delta-function potential barrier (see, 
e.g. [ pT|JT(i| ). The matching results in following linear 
system determining a, b, c and d 
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(9) 



Though complicated, the exact scattering probabilities 
can now be found and numerical results for different val- 
ues of the effective chemical potential have been given by 
Sipr and Gyorffy JT^ ]. However, the calculations may be 
simplified significantly for materials with a high Fermi 
energy compared to the temperatures or bias voltages of 
interest. In this limit transport only takes place near the 
Fermi level. 

We choose polar coordinates and allow for a Fermi ve- 
locity mismatch by considering a wave vector on the nor- 
mal side given by k <N) = fcp N) (sin 9 cos 4>, sin 9 sin (f>, cos 9) 
and a wave vector on the superconducting side with 
|k< s >|=fc F S) . 

The boundary conditions are satisfied only if fc<, N) = 
k£ } and fc^ N) = fc^ s) , as dictated by the transla- 
tional invariance along the interface. This means 
that kf ] = kp^i-rl sin 2 e, where the Fermi mo- 
mentum ratio is given by 77- = kp /kp . The 
wave vectors on the normal conducting side can now 
be written as q ± = fcp N) y/ cos 2 9 ± E//i <N) and on 
the superconducting side we similarly get fc ± = 



fc£° ^/(l - r\ sin 2 9) ± y/W- Af,/^ s \ This way of in- 
cluding the angle dependence is in accordance with Refs. 
P,^| JTllJl^ , [i^ , [i8|| but differs from the approach of Chaud- 
huri and Bagwell j?J in which the angle corrections to 
the ID expressions for q and k are approximated by 
cos# projection-factors. At low temperatures £7/i <N) ~ 
^/b 2 -Ao//i (S) <C 1, and therefore, we apply the Andreev 
approximation: k + = k~ = k ( z s > and q + = q~ = k^ . 
For values of /i (N,s) relevant for normal metals and low- 
temperature superconductors this approximation is valid 
for angles of incidence 9~ir/2. Semiconductors have 
much lower Fermi energies as compared to those of nor- 
mal metals, but, even for T/Tp ~ 1/50 the approxima- 
tion is reasonable. For angles in the vicinity of 9 ~ it/ 2 
the approximation becomes inaccurate. However, quasi- 
particles with vanishing perpendicular momentum do not 
contribute significantly to the perpendicular current and 
their effect in the IV curves and related quantities will 
thus not be important for semiconductors either. We 
have in fact checked this for the excess current by nu- 
merically solving Eq. (16) and found less than a half 
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TABLE I. Scattering probabilities at NS interface as a function of the normalized excitation energy E = E/Aq. 



A(E,l 
B(E, ( 

C(E,( 

D(E,( 



\e\ < i, e < e c 

1 

B2 + (l-B2)( 1 + 2z 2 ff(fl) ^ 

i-A(E,e) 





Igj > i, g < e c 

i 

[tB|+yi^T(l+2Z o 2 ff (9))J T 
4( J E 2 -l)(Z c 4 ff (9)+z; ff (9))^ 

2 V-B 2 -l(|-B| + Vs 2 -l) (l+Z 2 ff (0)) 
[|S|+V-B 2 -l(l+2Z c 2 ff (9))] 2 



2VB 2 -l(|B|-V / -B 2 -l)z c 2 ff (9) 
|B| + V£ ;2 -l(l+2Z 2 ff (9)) 



percent deviations for the GaAs 2DEG considered in 
FIG. 2. With these approximations the amplitudes be- 
come 

U V 
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v cos y 



Uq 



' 1+Vr 



cos I 



d = 



iv (r-^/r 



• l-rr> 



COS 



(10) 

(11) 

(12) 
(13) 



where we define 7 = + (uq — Vq) Z% s (0), and where 



z cS (6) = Wr(fl) 



cos 6* 



(£(gK - 1) 2 

4r(0)r„ 



(14) 



is an effective barrier strength, r v = Vp /Vp is the Fermi 
velocity ratio and T(9) = cos 9/y/i-rl sin 2 e. The dimen- 
sionless barrier strength Z = HJh^/v^K^ was intro- 
duced by Blonder and Tinkham ||. 

In order to obtain the scattering probabilities A, B, C 
and D we use the conservation of the probability-current- 
density, Eq. (||). For the z-direction this yields 



J 



\a^+\b^+e(\E\-A : 

A B ^ 



\ u 0\ 



Tr 



| c | 2 + e(|i?|-A ) 



Tr 



(15) 



It turns out that the scattering probabilities of the 
BTK model can still be applied provided that the dimen- 
sionless barrier strength is replaced by the introduced 
effective barrier strength given in Eq. (jl4|). For perpen- 
dicular incidence (6 — 0) this result agrees with the BTK 
result (r v = 1) M and the Blonder-Tinkham result @] 
which includes the possibility of a Fermi velocity mis- 
match. For a general angle of incidence and matching 
Fermi velocities and Fermi momenta (r v = = 1) the 
result reduces to that obtained by Kupka || . 

As mentioned, the wave vectors of the transmitted 
waves have the form kf } = kp \/i-rl sin 2 e. The square 
root defines a critical angel of incidence 8 Cl above which 
the solutions are evanescent and below which we have 
propagating waves, i.e. 6 C = arcsin(l/Vfc) for > 1. 
The physical reason for the critical angle is that the par- 
allel momentum exceeds the Fermi momentum of the su- 
perconductor and thus momentum can not be conserved. 
For rfc < 1 there is no critical angle due to the parallel 
momentum not being conserved. Going beyond the An- 
dreev approximation introduces another energy depen- 



dent critical angle 8 C = arcsin y/l - E/n^ caused by 
the wave vector q~ of the Andreev scattering state being 
imaginary JItJ . However, the later critical angle has little 
consequences for our results for the same reasons as when 
we discussed the validity of the Andreev approximation. 
In Sec. |V|, we suggest how the angle dependence may be 
probed. 

The directions of the reflected and transmitted waves 
can be obtained by considering the probability-current- 
density and the result is sketched in Fig. |l|. The angle of 
reflection 9 r coincides with the angle of incidence 6 and 
the angle of transmission is given by 

sin 6*4 = rfc sin# (16) 

in analogy with Snell's law in optics as it was also found 
by Kupriyanov |l(| . The general results for the scattering 
probabilities including the possibility of a Fermi velocity 
mismatch are summarized in Table |[ We conclude that 
the scattering probabilities of the BTK model still ap- 
ply, provided that the dimensionless barrier strength is 
replaced by an angle dependent effective barrier strength. 
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As the angle of incidence is increased, we observe an in- 
creasing effective barrier strength and therefore Andreev 
reflection is suppressed when the parallel momentum be- 
comes much larger than the perpendicular momentum 
pL In the same way, normal reflection increases when 
the parallel momentum increases. Nevertheless, we still 
have unit probability for Andreev reflection at the gap 
edge for all angles of incidence (9 < 9 C ) and for all bar- 
rier strengths. 





normal reflection 



Andreev reflection 



IV. CURRENT- VOLTAGE CHARACTERISTIC, 
EXCESS CURRENT AND DIFFERENTIAL 
CONDUCTANCE 



We calculate the current on the normal side of the in- 
terface where the current is carried only by single quasi- 
particles and no supercurrent. The current density in the 
z-direction is given by 



d d k 

(2n) d 



eve 2 f'(k) 



(17) 



where d — 1,2,3 is the dimension of the electron-gas 
and / <N) (k) is the non-equilibrium distribution function 
on the normal side of the interface. This approach ne- 
glects coherent effects of the propagation of electron-like 
and hole-like quasiparticles in the normal region and it 
applies to NS interfaces with a ballistic normal region 
and/or NS interfaces where the length of the normal re- 
gion is large on the scale of the phase coherence length. 
The integration is performed using polar coordinates ap- 
propriate for a ID electron gas, a two-dimensional (2D) 
electron gas, and a 3D electron gas, respectively. The ID 
case corresponds to the BTK model. 







transmission without branch-crossing 



transmission with branch-crossing 



FIG. 1. Directions of the transmitted and reflected waves 
in real space. The full dots (•) represent quasiparticles of 
predominantly electron-like character and the open dots (o) 
represent quasiparticles of predominantly hole-like character. 

In general the non-equilibrium distribution function 
can be found from a suitable transport equation, e.g. 
Boltzmann equation. Instead of taking this path, we 
follow BTK and assume that all quasiparticles incident 
from the reservoir are distributed in accordance with 
the Fermi-Dirac equilibrium distribution function with 
a shift in energy due to the applied voltage. When cur- 
rent flows, the reservoir is not in true equilibrium. How- 
ever, the voltage drop across the normal region can be ac- 
counted for by an Ohmic series resistance. We calculate 
f {N) {E,9) = e(TT/2-9)f^(E 7 9) + Q(9-TT/2)f^(E,9) 
by considering the two sub-populations separately. If wc 
take the chemical potential of the superconductor as ref- 
erence, we get f^ } {E,9,V) = f (E - eV) for the sub- 
population of quasiparticles with a positive momentum 
in the z-direction. The sub-population of quasiparticles 
with a negative momentum is 



J 



(E, it - 9, V) = A(-E, 9) [1 - f™(-E, 9, V)} + B(E, 9)fZ } (E, 9, V) 
+C(E, 0)f*> (E, 9 S , V) + D(E, 6)f™(E, 9 S , V). 



(18) 



Here the first term represents Andreev reflection of time- 
reversed quasiparticles, the second term represents nor- 
mal reflection, and the last two terms represent trans- 
mission of quasiparticles from the superconductor where 
f^(E,9 s ,V) = f (E) and 9 S = arcsin (r k sin0). Using 
the translational invariance along the interface, the sum 



I 

rule 1 = A + B + C + D, and the symmetries with respect 
to energy yields the normalized current 



Ao 



dET(E) f (E-eV/Ao)-f (E) , (19) 
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FIG. 2. Left panel: Normalized excess current cRnIcxc/ 'Ao(T) as a function of the dimensionless barrier strength Z. Right 
panel: Normalized differential conductance at zero bias g(V — 0) as a function of the dimensionless barrier strength Z. The thin 
lines, dashed lines and thick lines correspond to the ID BTK model, the 2D and 3D calculations, respectively. The temperature 
dependence of the excess current is entirely contained within the normalization through Ao(T) and the differential conductance 
is plotted at T — OK. The excess current and the zero bias differential conductance are shown for (a) parameters appropriate 
for a 2D GaAs-Al interface with a Fermi temperature Tf — 100 K of GaAs, (b) parameters appropriate for the 3D GaAs-Al 
interface studied by Taboryski et al. |l9| and (c) matching Fermi velocities and Fermi momenta. For the shown 2D result in 
(a) and 3D result in (b) the overestimation of Z for a perfect interface is 0.68 and 0.45, respectively. 
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where E = E/Aq is the normalized excitation energy, 
i?N is the normal state resistance and 



T 1D {E) 
T 2D (E) 

T 3D (E) 



1 + A(E,0)-B(E,0) 
Pna( 



= de 



2tt 



tt/2 

A9- 



1 + A(E,9) - B(E,( 

1 + A(E,9) — B(E, i 



(20) 
(21) 

,(22) 



are effective transmission coefficients for electrical cur- 
rent. The current for two and three-dimension systems 
have the same qualitative form as in the ID BTK model. 
However, quantitative changes are seen in the excess cur- 
rent and the differential conductance. 

The effective transmission coefficients are in general 
larger than the corresponding normal state transmission 
coefficients and this effect gives rise to a voltage depen- 
dent excess current compared to the normal state, where 
often the high-voltage limit 



eV»A 



I(V) 



lim I(V) 

An^O 



(23) 



is of interest from an experimental point of view. 

We have shown above, that a large parallel momen- 
tum suppresses the Andreev reflection probability and 
thus we expect to see a lower excess current in the three 
or two-dimensional limit as compared to the case of per- 
pendicular incidence. This is seen in the left panel of 
Fig. ^. For perfect [Z = 0) 2D and 3D interfaces with 
non-matching Fermi velocities and Fermi momenta, the 
ID BTK model overestimates the barrier strength Z sig- 
nificantly. For the shown 2D and 3D results the over- 
estimation of Z for a perfect interface is 0.68 and 0.45, 
respectively. 

At low temperatures the normalized differential con- 
ductance, g = Gns/Gnn, is given by 



g(V) = T(eV/A )/n 



(24) 



where Tn is the effective transmission probability when 
the superconductor is in the normal state. In the right 
panel of Fig. |^, results at zero bias are shown. Similarly 
to the excess current, we find a lower zero-bias conduc- 
tance with raising dimensionality as compared to the case 
of transport with perpendicular incidence. 

As an application of the present results, we now con- 
sider recent experiments by Taboryski et al. [p^[ who re- 
ported on Andreev reflections at interfaces between GaAs 
(3DEG) and superconducting Al films. The material pa- 
rameters are r v ~ 0.44 and ru — 0.03 and from the 
excess current of the ID BTK model, Taboryski et al. 
deduce the dimensionless barrier strength Zfi t to fall in 
the range from 0.7 to 0.9. Comparing with (b) in the 
left panel of Fig. || we find the barrier strength to fall 
in the range from 0.5 to 0.7. For GaAs, the energy de- 
pendence of the effective mass due to non-parabolicity 



is negligible within ~ 50Aai of the Fermi level ]2(J] for 
both the cases considered in FIG. 2. Since the zero-bias 
conductance is a Fermi-surface property no restrictions 
have been made by neglecting the energy dependence of 
the effective mass. For the high voltage limit of the 
excess current, the corrections due to a non-parabolic 
conduction band are small. For InAs, as considered by 
Schiissler and Kummel |p"l| the energy dependence due 
to non-parabolicity is more pronounced. 



V. SUGGESTED EXPERIMENT 

Benistant et al. |2^| have studied the angle dependence 
of Andreev scattering at Ag-Pb interfaces experimentally 
by using a magnetic focusing technique. The quasiparti- 
cles are injected to a very clean (ballistic) 3D Ag crystal 
through a point contact and the angle of incidence at the 
NS interface is controlled by a weak magnetic field. We 
suggest a variant based on an interface between a bal- 
listic two-dimensional electron-gas (2DEG) and a super- 
conductor and the technological opportunity of defining 
the angle of incidence geometrically. By applying gates 
on top of the 2DEG, it is possible to control the angle of 
incidence as sketched in Fig. [| 

Gale I Gate 2 Gate 3 



Superconductor 



FIG. 3. Interface between a ballistic 2DEG and a super- 
conductor. The three gates may individually be applied a 
gate-voltage and in this way it is possible to have perpen- 
dicular incidence or a finite angle of incidence defined by for 
example gate 2 and gate 3. 

If gate 1 and gate 3 are both negatively biased whereas 
gate 2 is turned off then quasiparticles have perpendicu- 
lar incidence. However, biasing gate 2 and gate 3, while 
gate 1 is turned off, a finite angle of incidence can be 
achieved. 

The angle dependence can be studied by measur- 
ing the differential conductance. At low tempera- 
tures the normalized differential conductance, g{eV, 9) = 
G NS (eV, 9)/G NN (eV, 6), is given by 



g(eV,9) = 



1 + T(9)A(eV, 9) - T(-0)B(eV, 9) 
\~T{-9)B{9) 



(25) 
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where B 



z? s /(i 



Z^a) is the normal reflection prob- 



ability when the superconductor is in the normal state. 
The transmission for quasiparticles leaving the interface 
through the gates, T(6), is peaked around 9 = 9o, with a 
width depending on the geometry. The possibility of mul- 
tiple Andreev reflections can be neglected if the phase- 
relaxation length is less than four times the distance be- 
tween the interface and the gate nearest to the interface. 

The experimental curves for g(eV, 0) and g(eV, 9$) may 
be fitted to Eq. (25) with the transmission, T(9q), and 
the barrier strength, Z, as fitting-parameters. Further 
information on Z and T(#o) m ay be obtained from the 
normal state conductance. When the superconductor is 
in the normal state the Landauer formula yields the con- 
ductance 



Gnn(#o) = ^-MT(9 ) 
h 



l-Bi 



(26) 



where M is the number of modes. Thus, agreement with 
the angle dependence is found if fits of the experimental 
curves for g(eV, 0) and g(eV, 9q) can be obtained using 
the same Z-value. 



VI. DISCUSSION AND CONCLUSION 

The angle dependence of scattering at NS interfaces 
is of important consequence when the NS interface has 
two or three-dimensional nature. The parallel degrees of 
freedom also have important consequences in supercon- 
ducting mesoscopic transport. If the current is carried 
by more than a single mode different modes represent 
different momenta parallel to the interface, and thus the 
scattering amplitudes depend on the mode index. More 
details may be found in a review of scattering theory in 
mesoscopic NS structures by Beenakker | Q . 

We have investigated the angle dependence of scatter- 
ing of quasiparticles at NS interfaces using the framework 
of Bogoliubov-de Gennes. As a main result the scatter- 
ing probabilities of the BTK model may still be applied 
provided that the scattering strength is replaced by an 
effective angle dependent barrier strength. This modi- 
fied effective scattering parameter agrees with previous 
calculations of BTK 0, Blonder and Tinkham M and 
Kupka Q. One of the consequences is that the Andreev 
reflection is suppressed for large angles of incidence and 
the normal reflection is increased towards unity. In the 
presence of a Fermi momentum mismatch, we find the 
angles of incidence and transmission to be related and 
in analogy with Snell's law, we find that above a certain 
critical angle of incidence we only have normal reflection. 

Furthermore, the results of the angle dependence have 
been applied to a NS interfaces with one, two and three- 
dimensional nature where we find that the ID BTK 
model overestimates the barrier strength. Calculations 
show that for certain material parameters and clean in- 
terfaces the corrections may be significant. However, the 



over-all qualitative predictions of the ID BTK model are 
found to agree with the new calculations. 
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